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ABSTRACT 

Two-dimensional numerical simulations of the effect of background turbulence on 2D 
resistive magnetic reconnection are presented. For sufficiently small values of the re- 
sistivity (77) and moderate values of the turbulent power (e), the reconnection rate is 
found to have a much weaker dependence on r\ than the Sweet-Parker scaling of 77 1 / 2 
Qh and is even consistent with an 77-independent value. For a given value of 77, the depen- 

q dence of the reconnection rate on the turbulent power exhibits a critical threshold in 

e above which the reconnection rate is significantly enhanced. 
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1 INTRODUCTION 

Magnetic reconnection is a ubiquitous plasma process 
whereby oppositely directed magnetic field lines break 
and rejoin in a different topological configuration — see, 
e.g., Biskamp ( 2000| ) for a review. It is probably the main 
mechanism behind many spectacular space- and astrophys- 
ical phenomena, such as sub-storms in the Earth's mag- 



netosphere ([Dungey 1961) and solar/stellar (Yokoyama et 
al. 2001) and accretion disk flares (Goodson et al. 1999). 
Reconnection has also been suggested as a possible mech- 
anism behind a number of high-energy astrophysics pro- 
cesses, such as magnet ar giant flares ( |Lyutikov|2003 ) , 7-ray 
bursts ( |Giannios &; Spruit||2006| ) and rapid TeV flares in 
blazars ( Giannios et al.|2009| . It is believed to lead to heat- 
ing and non-thermal particle acceleration in astrophysical 
coronae ( Drake et al.| |2006) and to play an important role 
in magnetohydrodynamic (MHD) turbulence and large-scale 
dynamos ( |Brandenburg Sc Subramanian||2005|. It is also a 
key element in sawtooth crashes in tokamaks ( |Hastie|1998 ). 

Early attempts to understand reconnection in terms of 
the simplest possible description of the plasma — single- fluid 
resistive MHD — led to the Sweet-Parker (SP) model ( | Sweet 
1958; Parker|1 957), which predicts a reconnection rate pro- 
port ionaltoS^^^, where S — LVa /rj is the Lundquist num- 
ber, L is the system size, Va is the Alfven speed and 77 is 
the plasma resistivity. In many astrophysical environments, 
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S ^> 1 [e.g., S ~ 10 14 in solar flares, S ~ 10 18 in the in- 
terstellar medium (ISM)], resulting in SP reconnection rates 
which are orders of magnitude slower than that observed in 
the above mentioned phenomena. The main quest in recon- 
nection research has thus been to identify the key physical 
mechanisms that are missing from the SP theory and that 
can explain these fast rates. Following much numerical work, 
it is now thought that in low- density, collisionless plasmas, 
nonclassical effects [e.g., the Hall effect ( Biskamp |2000| or 



anomalous resistivity ( Malyshkin|2005 )] can indeed give rise 
to fast, Petschek-like iPetschek 1964), reconnection. Fast, 



collisionless reconnection can only take place when the re- 
sistive width of the reconnection layer, 5sp ~ L/S 1 ^ 2 (Sweet 
1958 Parker 1957) is small compared to the relevant ki- 



netic scale [which is frequently the ion collisionless skin- 



depth (Cassak et al 
|2007| 



2005 Yamada et al. 2006 Uzdensky 



di 



— c/ujpi, or else d e — c/u pe in the case of pair 
plasmas]. In many astrophysical situations, e.g., the solar 
chromosphere, the ISM, inside stars and accretion disks, 
and in the high-energy-density environments in central en- 
gines of gamma-ray bursts (GRB) and core-collapse super- 
novae ( |Uzdensky &; MacFadyen|2006| ) , the density is so high 
that the above condition is not satisfied; for example, for 
the warm ionized ISM, 5sp/di ~ 10 3 [and even larger for 
diffuse clouds or molecular clouds (Zweibel 1989)]; in the 
magnetosphere of the accretion disk in the collapsar sce- 



nario of long GRB central engines ( Uzdensky & MacFadyen 
2006) Ssp/de ~ 10 7 (here d e may be more important than 



because, under some circumstances, this is expected to 
be mostly pair plasma). Therefore, the reconnection layer 



2 N. F. Loureiro et al 



is collisional and resistive MHD should apply. The prevail- 
ing opinion is, however, that in the resistive-MHD case the 
Petschek mechanism fails (Biskamp 1986; Uzdens ky &; Kul^| 
|srud|2000| , and reconnection is negligibly slow, possibly de- 
scribed by the SP theory. Can fast (i.e., 77— independent) 
reconnection happen in these environments? 

Missing from the SP picture is that in Nature, most, 
if not all, environments where reconnection takes place are 
likely to be turbulent [e.g. (Retino 2007)]. Can background 
turbulence significantly accelerate reconnection? The aim of 
this Letter is to address this question. 

Turbulent reconnection has been studied previously, 



1986 


|1988| |Lazarian & Vishniac| [1999 ; Kim & Diamond 


2001 


and numerically (|Matthaeus & Lamkin 1986 


Fan et 


al. |2004| |Smith et al.| 


2004; |Watson et al.||2007| |] 


^apenta 


2008 


Kowal et al.|2009 


) . Given the complexity of the prob- 



lem, however, a rigorous theory of turbulent reconnection is 
not yet available and numerical studies have been limited 
by resolution constraints. Thus, the role of turbulence in 
reconnection remains controversial. A pioneering numerical 
study with decaying 2D turbulence by Matthaeus &; Lamkin| 



(1986) suggested that a finite level of broadband MHD fluc- 



tuations can enhance reconnection [see also Matthaeus & 



Montgomery (1986)], but the very limited computing ca- 
pabilities employed precluded a clear conclusion about the 
asymptotic regime n — > 0. Subsequently, in a highly influ- 



ential paper, Lazarian & Vishniac (hereafter, LV) ( |Lazar- 
ian &; Vishniac|1999 ) suggested that turbulence can greatly 



accelerate reconnection by enabling multiple reconnection 
sites in the current sheet. The LV picture is an essentially 
three dimensional (3D) process, as the multiple reconnec- 
tions mechanism they envision is topologically prohibited in 
two dimensions (2D). Recent 3D numerical simulations by 
Kowal et al. ( 2009| seem to support the LV model; however, 
due to computational constraints, only moderate scale sep- 
arations and fairly strong turbulence (but still sub-Alfvenic 
compared to the reconnecting field) could be probed. 

In 2D, the LV mechanism should not work and the ex- 
pectation has been that topological constraints should pre- 
vent significant acceleration of reconnection by turbulence. 
In this Letter, we present evidence contrary to this belief: 
our numerical results show that background turbulence can 
have a dramatic effect on 2D magnetic reconnection, yielding 
reconnection rates that exhibit a much weaker dependence 
on S than the SP scaling of S -1 ^ 2 and may even asymp- 
tote to an ^-independent value as S — > 00. These results 
call for a further theoretical effort to understand magnetic 
reconnection in the presence of turbulence. 



2 NUMERICAL SETUP 

Our aim is to investigate how relatively weak MHD turbu- 
lence affects reconnection in 2D. The incompressible MHD 
equations, with an external forcing term F, are: 

dtu + u • Vtt = B • VB - V (p + B 2 /2) + v\J 2 u + F, (1) 



d t B + u • VB = B Vu + nV 2 B, 



(2) 



where u is the velocity, B is the magnetic field, and p is 
the pressure. The pressure gradient is determined by the in- 
compressibility condition V • u — (which sets the plasma 



density p — const), so an energy equation is not required. 
The resistivity of the plasma is denoted by 77 and the viscos- 
ity by v. We have normalized p and B by p and (47rp) 1//2 , 
respectively. Expressing the 2D [(x, y)] incompressible ve- 
locity field (V • u = 0) in terms of the stream function (f) 
and the solenoidal magnetic field (V • B = 0) in terms of 
the flux function ijj, u — (—d y (j),d x 4>), B = {—d y i^,d x ^), 
we can rewrite the above equations in terms of <j) and ip as 
follows ( |Strauss||1976| ): 



6> t V 2 0+{0, V 2 0} 

6>W + {<^} 



{^vivj + ^vV + Z, (3) 

77V 2 (i/> - ^ eq ) , (4) 



where the Poisson brackets are denoted by {0, ip} = 
dx^dyij) — d y (j)d x ^. These are the equations that we solve 
with our numerical code. The dimensions of length are 
set by the box size (L x ,L y ) = (27r, 2.187r) and the di- 
mensions of magnetic field are set by the equilibrium con- 
figuration: we use a tearing-mode unstable equilibrium, 
V^q = i/jo I cosh 2 [x) , where ipo = 3>/3/4, so that the 
maximum value of the equilibrium magnetic field, Bo — 
(d'0eq/dx) m ax = 1. The background density is normalized 
to 1 and the dimensions of velocity are given by setting 
the Alfven speed Va = Bq = 1. Thus, time is measured 
in Alfvenic units, L x /(2tyBo). Diffusion of the equilibrium 
magnetic field is prevented by the addition of an external 
electric field to the RHS of Eq. ^ — the term -ryV^eq- 
The externally imposed (turbulent) forcing / = z • (V X F) 
(random, white noise in time), is characterized by two pa- 
rameters: e = (u • F), the power input per unit area and kf, 
the forcing wave number. Note that the forcing is applied 
only to the momentum equation, and does not break the 
frozen- flux constraint. The equations are solved in a dou- 
bly periodic box using a pseudo-spectral, symplectic algo- 



rithm (Loureiro &; Hammett||2008| (which can be used in 



either semi-implicit or explicit timestepping mode: the for- 
mer is employed in the laminar simulations; the latter in the 
turbulent ones). Our ultimate goal is to characterize the ba- 
sic behaviour of the effective reconnection rate E e f? (defined 
below) in the 4D parameter space: E e s(rj, e,kf). However, 
in the present study we focus on only two parameters, r\ — v 
(i.e., the magnetic Prandtl number Pm = is/rj = 1) and e. 

We start our simulations without turbulence follow- 
ing a standard tearing mode evolution. After the lin- 
ear ( |Furth et al.|1963|), |Rutherford| ( fl973| , and X-point col- 
lapse stages ( Loureiro et al.||2005 ), a thin SP current layer 
forms between two large magnetic islands. For our parame- 
ters, this happens when the value of the reconnected flux at 
the X-point is ipx ~ 0.4, whereupon the length of the cur- 
rent sheet is L w 1, so we define S = LVa/t] = 1/rj. Since 
our intention is to focus solely on the effect of turbulence 
on resistive reconnection, we use the laminar SP configu- 
ration at ipx ~ 0.4 as a starting point for our turbulent 
runs and switch on the driving term / in Eq. (|3| at this 
stage. Furthermore, when ipx ~ 0.9, the SP stage gives way 



to saturation (Loureiro et al. 2005). Thus, we restrict our 



analysis of the numerical data to the time interval where 
0-4 < < 0.9. 

Our runs are characterized by the following parameters. 
We scan the range in e from e = 3x 10 -4 to e = 0.1, and in n 
from n = 10 -3 to n = 6.5 x 10 -5 . The forcing wave number is 
kf = 10, giving turbulent motions with characteristic scales 
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Figure 1. Magnetic reconnection in the presence of background 
turbulence. Contour plot of current density j z = V 2 ^ at t = 735 
from the run with parameters 77 = 6.5 x 10 — 5 , kf = 10 and 
e = 0.01 (cf. blue dotted line of Fig. [3}. Run resolution is 8192 2 . 

a few times smaller than the current sheet length (If ~ 0.3). 
For 77 = 6.5 x 10 -5 , the width of the laminar current sheet 
is Scs ~ 0.014, so we have a reasonably good scale sepa- 
ration in these simulations. Resolutions up to 8192 2 were 
used. Such resolutions would not have been possible in 3D. 
Convergence tests were performed on selected runs to ensure 
that increasing the resolution did not change the results. 




Figure 2. Contour plots of the current density at different times 
for the runs with 77 = 6.5 x 10 -5 , kf = 10 and e = 0.01 (top row, 
cf. blue dotted line of Fig. [3} and e = 0.1 (bottom row, cf. red full 
line of Fig. [3}. Resolution is 8192 2 . Selected contour lines of ip 
are overplotted. Snapshots zoom in on the current sheet region. 



3 RESULTS 

A typical snapshot of these simulations is shown in Fig. [I] 
The large-scale reconnecting configuration is manifest; the 
turbulent motions, which are also clearly visible, are not suf- 
ficiently strong to destroy itQ so the reconnection problem 
is still well posed. In the presence of turbulence, the current 
sheet is no longer straight as in the laminar case, but now 
wiggles about its original (i.e., laminar) location (x = 0): 
see Fig. [2] We find that secondary islands (plasmoids) are 
constantly present, being continuously generated and ex- 
pelled from the sheet (whereas the laminar simulations at 
most display one plasmoid). The size of these plasmoids in- 
creases with e, as does the turbulent distortion of the sheet. 

We would like to measure the reconnection rate, E e f£, 
in the presence of turbulence. In the laminar case, E e R = 
dipx/dt. However, in turbulent simulations, the wandering 
of the current sheet makes the reconnection rate difficult to 
define. To quantify the effective reconnection rate, smoothed 
over the rapid turbulent fluctuations, we employ the fol- 
lowing procedure. Consider the representative timetraces of 
ip(x = 0, y — 0) of Fig. [3] We see that they can be described 



1 There are regions with qualitatively different turbulence. This 
is because, whereas our turbulent forcing is statistically homo- 
geneous, the background magnetic and velocity fields associated 
with the large-scale reconnecting configuration are not. 
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Figure 3. Selected time traces of the magnetic flux at the location 
of the laminar X-point (x = 0, y = 0), for 77 = 6.5 x 10 -5 , kf = 10 
and different values of e. 



by mostly upward fluctuations on top of a monotonically in- 
creasing baselin^] We have checked that the minima of these 
curves correspond to maxima of the current. The reason is 
that, because of the turbulence, the actual very thin current 



2 Occasional downward fluctuations are due to plasmoids cross- 
ing the (x = 0, y = 0) location. 
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Figure 4. Plot of the effective reconnection rate (electric field) as 
a function of the Lundquist number S = 1/rj for different values 
of e and kf = 10. Dashed line shows S~ 1//2 scaling. Data for 
e = 0.1 and S ^ 3 x 10 3 is not shown because the magnitude 
of the fluctuations and the relatively low values of S prevent a 
statistically accurate computation of E e &. 



sheet wanders around the original (laminar) X-point loca- 
tion — see Fig. [2] As it deviates from this location, the value 
of ipx increases, since the laminar ip profile has a minimum 
at the X-point. When the current sheet crosses the location 
of the laminar X-point, the measured current is maximum, 
and ip is again minimum. This suggests that the slope of a 
linear fit to the local minima of these timetraces is an ac- 
curate diagnostic of the reconnection rate and that is what 
we use. Error bars are estimated by reducing the interval of 
values of ip where the minima are selected by 0.05 at both 
ends. Although the turbulence is switched on at ipx ~ 0.4, 
we only apply this diagnostic for ip > 0.5 so that the turbu- 
lence can reach a saturated stage before we assess its effect 
on the reconnection rate. We emphasize that this is not the 
only possible diagnostic of the reconnection rate and indeed 
we have tried several different ones, with similar results. 

The resulting dependence of the effective turbulent re- 
connection rate E e s on S = 1/rj is shown in Fig. [4] for a 
range of values of the injected power e. For comparison, we 
also plot the reconnection rate obtained in the laminar runs, 
Esp, which yields very good agreement with the SP scaling 
of S -1 / 2 (dashed line). The immediate observation, and the 
main result of this Letter, is that, at sufficiently large values 
of S, all sequences with e > 0.003 exhibit a dependence on 
S that is significantly shallower than the SP dependence of 
S -1 / 2 and may even consistent with a S -independent value 
as S — > 00. Also remarkable is that, common to all finite- e 
sequences, there is a clear transition that does not seem to 
depend on e, always taking place around S ~ 2 x 10 3 . Specif- 
ically, at values of S below this, the reconnection rate in the 
turbulent runs is enhanced over the SP value, but only mod- 
estly, and still scales as ~ S -1 ^ 2 . Above this critical value 
turbulence has a much stronger effect on reconnection. 

Shown in Fig. [5] is the dependence of the reconnection 
rate enhancement factor E e ^/Esp on the injected power e. 
There is a distinct transition at e ~ 0.003. For smaller e, 
the dependence of the reconnection rate on e is very weak, 
roughly consistent E e ^/Esp ~ 1; for e > 0.003, the enhance- 
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Figure 5. Plot of the reconnection rate enhancement factor as 
a function of the injected (turbulent) power, e, for the smallest 
values of 77 and kf = 10. 



ment factor increases with e in a way which is consistent both 
with a power law E e ^/Esv ~ e a , where a ~ 0.15 — 0.25, and 
with a logarithmic dependence. Note that the power law ex- 
ponent is significantly shallower than the LV prediction of 
e 1 / 2 for the 3D case. 



4 DISCUSSION AND CONCLUSION 

In summary, we have performed a detailed numerical study 
of the effect of background turbulence on 2D incompressible 
resistive magnetic reconnection. Our key finding is that, at 
sufficiently large values of the Lundquist number and mod- 
erate values of the injected (turbulent) power, the scaling of 
the reconnection rate with S is much weaker than the SP 
scaling of S -1 ^ 2 and is, in fact, consistent with asymptoti- 
cally fast reconnection (independent of S; see Fig.j4|. Other 
important conclusions are: (i) the transition from slow to 
fast reconnection happens at a fixed (i.e., independent of e) 
value of the Lundquist number, S w 2 x 10 3 and (ii) the 
enhancement of the reconnection rate over its laminar (SP) 
value exhibits a clear threshold in the injected turbulent 
power, e w 0.003, below which turbulence has little influ- 
ence on the reconnection rate and above which significant 
enhancement of the reconnection by turbulence starts — 
see Fig. |5] In our simulations, we find that the rms turbulent 
velocity u rms (averaged over the simulation box), is closely 
described by ?x rms « 2.44(e/A;/) 1/3 . Thus the Alfvenic Mach 
number associated with the transition to fast reconnection 
is Ma = Utuis/Va ~ 0.16 (note that the dependence of our 
results on kf has not yet been tested; this can influence the 
numerical values of both the S and the e transitions) . 

Our results imply that a novel mechanism of turbulent 
enhancement of reconnection exists which is operative al- 
ready in 2D. While our results support the basic claim by 
LV of turbulent enhancement of reconnection, we empha- 
size that, because they are 2D, they cannot be explained by 
any present theory, including theirs. One speculative pos- 
sibility is that the enhancement of the reconnection rate is 



related to the formation of multiple plasmoids (Loureiro et 



al. 2007 2009 Samtaney et al. 2009). This hypothesis is 
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strengthened by the fact that the transition from slow to 
fast reconnection (see Fig. |4| is independent of e and oc- 
curs at a value of S consistent with the critical Lundquist 
number for plasmoid generation ( |Loureiro et al.|2005 ). As S 
increases, turbulence might facilitate frequent plasmoid for- 
mation by the local, transient, enhancement of the tearing 
mode instability parameter A' ( Furth et al.| [l963) (i.e, the 
local enhancement of the aspect ratio of the current sheet), 
but whether plasmoids can indeed provide the enhancement 
of the reconnection rate that we observe is not yet known. 

In this study, the effect of viscosity has not been ad- 
dressed. Our modelling choice has been to keep the magnetic 
Prandtl number fixed and equal to 1. This is not character- 
istic of most astrophysical systems, where either Pm ^> 1 
(e.g., warm interstellar and intracluster medium, as well as 
some accretion disks) or Pm <C 1 (e.g., stars, planets and 
liquid-metal laboratory dynamos) (iBrandenburg &; Subra- 



2005). We have adopted Pm — 1 to avoid addi- 



tional complications and to present a clear demonstration 
of principle — that turbulence can accelerate reconnection 
in two dimensions. Changing Pm from 1 will result in two 
effects: i) different viscous cutoff scale; however, since our 
results appear to depend only on the outer scale features 
of the turbulence (i.e. e, and possibly also kf), we do not 
expect this effect to be significant; ii) different laminar solu- 
tion, with a broader current sheet and slower outflows (Park 
et al. || 1984 ). This might affect the instability of the current 
sheet to plasmoids and could potentially impact on the re- 
connection rate in the presence of turbulence. Investigations 
of how Prandtl number (and forcing scale) might affect our 
results are left for future work. 

A final caveat is that, in this work, we have neglected 
two- fluid effects, whose importance in reconnection is now 
widely recognized. There are many astrophysical environ- 
ments where these effects will doubtlessly play an important 
role (e.g. solar flares); a future publication will assess their 
importance in turbulent reconnection. 
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